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A  previous  theory  [M.  Ershov  et  al.,  Appl.  Phys.  Lett.  67, 3147  (1995)]  for  studying  the  distribution 
of  nonuniform  fields  in  multiple-quantum-well  photodetectors  under  an  ac  voltage  is  generalized  to 
include  nonadiabatic  space-charge-field  effects.  From  numerical  results  calculated  by  the 
generalized  theory,  it  is  found  that  field-domain  effects  are  only  important  at  high  temperatures  or 
high  voltages,  where  both  injection  and  sequential-tunneling  currents  are  expected  to  be  large.  On 
the  other  hand,  field-domain  effects  become  negligible  at  low  temperatures  and  low  voltages,  but 
nonadiabatic  effects  included  in  this  extended  theory  are  enhanced  for  small  sequential-tunneling 
currents.  The  time  duration  for  nonadiabatic  effects  is  determined  by  the  quantum  capacitance.  By 
using  the  generalized  theory,  a  differential  capacitance  is  calculated  for  a  non-steady  and  a 
negative  conduction  current  is  predicted  under  a  positive  voltage  in  this  case  due  to  charge 
accumulation  around  the  collecting  contact. 

[DOI:  10.1063/1.1594815] 


I.  INTRODUCTION 

Multiple-quantum- well  (MQW)  photodetectors  using  in¬ 
tersubband  transitions  have  attracted  a  lot  of  studies  over  the 
past  few  years.1  Transient  spectroscopy  allows  us  to  gain 
information  on  the  sequential-tunneling  processes  between 
quantum  wells  (QWs)  and  thus  QW  parameters,  including 
geometrical  and  QW  capacitances  at  the  same  time. 

When  a  voltage  is  applied  to  a  MQW  structure,  electrons 
respond  by  producing  a  sequential-tunneling  current.  The  in¬ 
teresting  thing  is  that  a  uniform-field  distribution  inside  the 
MQW  structure  is  not  stable  if  injection  and  sequential- 
tunneling  currents  are  quite  different  from  each  other.2  With 
the  existence  of  an  imbalance  between  injection  and 
sequential-tunneling  currents,  the  uniform  electric  field  is 
split  into  many  local  fields.  The  local  electric  fields  which 
are  constant  inside  barrier  layers  are  all  different  from  each 
other,  leading  to  a  distribution  of  nonuniform  electric  fields 
within  the  MQW  structure.  The  distribution  of  nonuniform 
fields  depends  on  the  applied  voltage,  and  is  accompanied  by 
a  nonuniform  distribution  of  electron  densities  in  different 
QWs.  When  the  voltage  varies  with  time,  the  distribution  of 
nonuniform  fields  evolves  with  time,  giving  rise  to  varying 
field  domains.  The  difference  between  two  local  electric 
fields  within  neighboring  barrier  layers  is  determined  by  the 
charge-density  deviation  from  the  equilibrium  value  inside 
the  sandwiched  QW  through  the  boundary  condition  contain¬ 
ing  the  quantum-well  capacitance. 

*  The  MQW  capacitance  contains  two  parts.  The  first  part 
is  the  geometrical  capacitance  which  depends  on  the  dielec¬ 
tric  constant  of  the  host  semiconductors.  The  second  part  is 
the  quantum  capacitance  which  is  attributed  to  the  density- 
of-states  of  the  two-dimensional  electron  gas  in  each  QW. 
The  geometrical  capacitance  is  associated  with  classical 
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charge  fluctuations2  after  an  ac  voltage  is  applied.  However, 
the  quantum  capacitance  is  related  to  the  quantum  charge 
fluctuations  resulting  from  nonadiabatic  space-charge-field 
effects.3,4  The  nonadiabatic  nature  of  these  effects  is  due  to 
the  fact  that  the  transient  currents  that  appear  depend  not 
only  on  the  time-dependent  electric  field,  but  also  on  its  time 
derivative. 

Most  of  the  previous  calculations  on  the  distribution  of 
field  domains  were  limited  either  to  dc  voltage5  or  to  ac 
voltage  in  steady  state.  Very  few  calculations  were  done  for 
ac  voltage  in  non-steady  state.4  In  steady  state,  charge  fluc¬ 
tuations  are  adiabatic,  namely  the  time  derivatives  of  the 
charge  densities  in  the  QWs  at  each  moment  are  zero.  The 
detected  conduction  current  equals  the  injection  current 
which  is  modified  by  the  local  field  at  the  emitter  barrier. 
However,  the  time  derivatives  of  the  charge  densities  are  no 
longer  zero  in  non-steady-state  conditions.  In  this  case, 
charge  might  be  dynamically  accumulated  or  depleted  in  the 
MQWs.  Adiabatic  charge  fluctuations  only  depend  on  the 
instantaneous  voltage  at  each  moment,  while  nonadiabatic 
charge  fluctuations  depend  on  both  the  voltage  and  its  time 
derivative. 

When  the  sequential-tunneling  current  is  low,  the  impu¬ 
rity  or  defect  channels  within  the  barrier  will  play  a  role.6 
However,  this  only  modifies  the  resistance  of  the  sample  for 
the  sequential  tunneling  of  electrons.  The  nonadiabatic  ef¬ 
fects  discussed  in  this  article  for  electron  tunneling  remain 
the  same.  For  the  same  reason,  the  variation  in  barrier  width 
and  alloy  composition  will  not  affect  nonadiabatic  effects 
discussed  in  this  article.  Moreover,  the  self-consistent  Har- 
tree  model  with  nonadiabatic  effects  on  non-steady-state 
electrons  under  an  ac  electric  field  in  multiple  quantum  wells 
has  proven  that  the  nonadiabatic  effect  discussed  in  this  work 
is  related  to  a  shaking  of  the  Fermi  level  in  the  nonequilib¬ 
rium  distribution  function  of  electrons  7 
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The  adiabatic  field-domain  theory2,5  has  been  proven  to 
be  a  great  success  to  explain  the  quantitative  difference  be¬ 
tween  the  experimental  and  uniform-field  theoretical  results 
at  T—ll  K  and  high  voltages.  On  the  other  hand,  the  nona- 
diabatic  tunneling  model4  under  a  uniform  electric  field  has 
also  been  demonstrated  successfully  in  explaining  the  experi¬ 
mentally  found  residual  current3  at  zero  ac  voltage  and  T 
=  40  K.  Therefore  it  is  quite  reasonable  to  expect  that  the 
combination  of  these  two  theories  will  provide  a  complete 
description  for  electron  sequential  tunneling  in  MQWs  at 
temperatures  either  above  77  K  or  below  40  K.  It  has  been 
found  that  the  nonadiabatic  behavior  (nonzero  residual  cur¬ 
rent)  observed  in  transient  tunneling  currents  as  a  function  of 
bias  voltage3  can  not  be  simply  explained  by  Levine’s  adia¬ 
batic  sequential-tunneling  model1  even  after  the  field-domain 
effect2  has  been  included.  Therefore  we  require  a  nonadia¬ 
batic  sequential-tunneling  model  with  the  inclusion  of  the 
field-domain  effect  to  explain  related  experimental  results  in 
a  wide  range  of  temperatures.  This  is  the  motivation  for  this 
work. 

In  this  article,  we  generalize  the  previous  theory2  for 
studying  the  distribution  of  nonuniform  fields  in  MQW  pho¬ 
todetectors  under  an  ac  voltage  by  including  nonadiabatic 
space-charge-field  effects.  We  find  from  numerical  results 
that  field-domain  effects  are  only  important  at  high  tempera¬ 
tures  or  high  voltages  when  both  injection  and  sequential- 
tunneling  currents  are  significant.  On  the  other  hand,  we  find 
that  nonadiabatic  effects  included  in  the  generalized  theory 
become  visible  at  low  temperatures  and  low  voltages  when 
field-domain  effects  are  negligible.  The  time  duration  for 
nonadiabatic  effects  is  found  to  depend  on  the  quantum  ca¬ 
pacitance,  while  the  classical  dielectric  displacement  current 
is  found  to  be  related  to  the  geometrical  capacitance.  More¬ 
over,  a  negative  conduction  current  is  predicted  under  a  posi¬ 
tive  voltage  in  non-steady  state.  This  article  only  discusses 
effects  of  classical  and  quantum  charge  fluctuations  on  the 
dark  current  in  multiple-quantum-well  photodetectors.  How¬ 
ever,  these  effects  on  photocurrent  can  be  discussed  in  a 
similar  way  if  we  assume  the  tunneling  current  in  this  case  to 
be  the  sum  of  the  dark  current  and  photocurrent. 

The  organization  of  the  article  is  as  follows.  In  Sec.  II, 
we  introduce  our  model  for  nonadiabatic  charge-density  fluc¬ 
tuations  and  the  distribution  of  field  domains  in  MQWs. 
Both  classical  and  quantum  charge  fluctuations  in  QWs  are 
formulated  and  studied.  Numerical  results  and  discussions 
are  given  in  Sec.  HE  for  effects  of  classical  and  quantum 
charge  fluctuations,  including  dependence  on  the  total  num¬ 
ber  of  quantum  wells  N,  temperature  T ,  capture  probability 
Pc ,  time  period  tp ,  and  amplitude  £0  of  an  applied  ac  elec¬ 
tric  field,  QW  electron  density  n2D ,  and  contact-layer  elec¬ 
tron  concentration  nc.  The  article  is  briefly  concluded  in 
Sec.  IV  with  some  remarks.  . 

II.  MODEL  AND  THEORY 

In  this  section,  we  first  study  the  effects  of  field  domains 
resulting  from  the  imbalance  between  injection  and 
sequential-tunneling  currents  and  show  that  these  effects  be¬ 
come  negligible  at  low  temperatures.  After  that,  we  general- 


FIG.  1.  The  processes  (a)-(d)  for  the  formation  of  field  domains  in  a 
multiple-quantum-well  (MQW)  sample  with  N  quantum  wells  (QWs).  Here, 
£k  for  fc=0,  1,  2,...,N  indicates  local  electric  fields  in  (V+ 1)  different 
barrier  layers,  and  pk  for  k~  1,  2 corresponds  to  different  charge  den¬ 
sities  inside  N  QWs.  is  the  chemical  potential  in  contact  layers,  Vb  is  the 
applied  voltage,  and  L t  is  the  total  length  of  MQWs. 

ize  the  adiabatic  field-domain  theory  by  including  nonadia¬ 
batic  space-charge  effects  resulting  from  quantum 
capacitance. 

A.  Classical  charge  fluctuations 

Electrons  in  QWs  are  confined  in  the  direction  perpen¬ 
dicular  to  the  wells,  while  electrons  in  heavily  doped  contact 
layers  are  free  in  all  three  directions.  Therefore  we  expect  the 
tunneling  of  electrons  from  a  contact  layer  to  a  QW  (3D-to- 
2D)  will  be  physically  different  from  that  of  one  QW  to 
another  (2D-to-2D). 

As  shown  in  Fig.  1,  we  see  that  the  distribution  of  uni¬ 
form  dc  electric  fields  in  (a)  with  £0  =  =  £^  =  ’  * ' =  £N  is  not 
stable  if  the  injection  current  flowing  from  the  left  contact 
layer  to  the  first  QW  and  the  sequential-tunneling  current 
flowing  from  the  first  QW  to  the  second  QW  are 
different.2,5,8  As  an  example,  we  assume  in  Fig.  1(b)  that  the 
injection  current  is  smaller  than  the  sequential-tunneling  cur¬ 
rent.  In  this  situation,  the  local  field  £0  in  the  first  barrier 
(emitter  barrier)  layer  has  to  be  increased  so  as  to  equalize 
these  two  currents.  As  a  result  of  £0>£i  =  E2=:'  "=£n>  we 
know  from  Maxwell  equations  that  the  charge  density  in  the 
first  QW  will  be  reduced  relative  to  the  others,5  i.e.,  p\<p2 
—  Pi  —  '"  —  Px .  Now,  let  us  further  compare  the  tunneling 
currents  flowing  from  the  first  QW  to  the  second  QW  and 
that  flowing  from  the  second  QW  to  the  third  QW.  We  realize 
from  Fig.  1(c)  that  the  tunneling  current  flowing  from  the 
first  QW  is  less  than  that  flowing  from  the  second  QW  since 
p  i  <  p2  for  £\  =  £2-  Therefore  £x  >  £2 = *  •  *  ==  £N  is  required  so 
as  to  equalize  these  two  sequential  tunneling  currents.  Con¬ 
sequently,  we  are  left  with  p j  <  p2 <  p3  —  •  •  ♦  =  pN .  This  pro¬ 
cess  will  continue  up  to  the  last  (Mb)  QW,  as  displayed  in 
Fig.  1(d),  until  the  initial  field  distribution  50>^1>£2>*** 
>£n  and  density  distribution  Pi<Pa<P3<' **<p^  are 
reached.  Here,  the  local  field  £j(t)  is  a  constant  in  the  jth. 
barrier  layer,  and  the  splitting  of  different  local  electric  fields 
constructs  the  field  domains  in  MQWs.  The  electric  fields  in 
different  barrier  regions  in  steady  state  will  be  redistributed 
if  the  total  tunneling  current  flowing  into  the  bottom  contact 
layer  is  different  from  the  injection  current  flowing  out  of  the 
top  contact  layer.  Therefore  the  above  field-adjustment  pro- 
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cess  will  be  repeated  again  and  again  until  the  total  tunneling 
current  to  the  bottom  contact  layer  is  equal  to  the  injection 
current  from  the  top  contact  layer,  and  then,  a  stable  distri¬ 
bution  of  local  fields  is  formed.  It  should  be  noted  that  the 
distribution  of  electric-field  domains  in  steady  state  depends 
on  the  value  of  voltage  applied  at  each  moment.  A  variation 
of  the  applied  voltage  with  time  results  in  varying  field  do¬ 
mains,  and  charge-density  oscillations  in  each  QW  as  well. 

The  injection  current  density  (3D-to-2D)  from  the  con¬ 
tact  layer  to  the  first  QW  is  calculated  to  be5 


4?I>o>^i(0,£oW] 


em*kBT' 


£b(0 


\+°°dE7[E,\£om 

Jo 


Xln 


1  +  exp[(fL0-E)/kBT] 


(1) 


where  t  is  the  time,  m  *  is  the  effective  mass  of  electrons,  T 
is  the  electron  (or  lattice)  temperature,  LB  is  the  thickness  of 
the  barrier,  l[E,\£(t)\]  is  the  quantum  transmission  of 
electrons6  with  kinetic  energy  E  through  a  barrier  biased  by 
an  electric  field  £(/),  and  the  chemical  potential  jnQ  in  the 
contact  layer  is  related  to  the  electron  concentration  nc  by 


-l 


For  low  T ,  electrons  in  QWs  can  only  populate  the  ground 
subband  with  quantized  energy  E0 .  Furthermore,  the 
sequential-tunneling  current  density  (2D-to-2D)  from  the  k\ h 
QW  to  the  neighboring  ( k+  l)th  QW  is  found  to  be1 


dE7[E+E0,\£k(t)\] 


X  -folE  -fik+i(t)  +  eLB\Sk(t )  |]}, 

W 

where  Lw  is  the  well  width,  k=  1,2 represents  the  in¬ 
dex  of  N  wells,  £k(t)  is  the  local  electric  field  in  the  Mi 
barrier  region,  k=NP 1  corresponds  to  the  bottom  contact 
layer,  fiN+i  =  ^Q,  fo(x)  =  [l^&^p(x/kBT)Yl  is  the  Fermi- 
Dirac  function,  and  the  chemical  potential  fik{t)  in  the  Mi 
QW  introduced  in  Eqs.  (1)  and  (2)  can  be  determined  by  the 
electron  density  nk(t)  in  the  Mi  QW  through 

[  irh2nk{t)\ 


Pk(t)  =  kBT\n\  exp 
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which  is  measured  from  E0  ■  In  Eq.  (2),  vk[8k (f)]  is  the  drift 
velocity  of  electrons  in  the  Mi  barrier  layer,  given  by1 


4^(0]= 


4(0 


jeHt)+e] 


with  the  saturation  velocity  and  saturation  field  £s ,  re¬ 
spectively. 

As  explained  in  Fig.  1,  the  distribution  of  field  domains 
changes  with  the  value  of  the  voltage  Vb(t)  applied  to  the 
sample.  When  Vb(t)  varies  with  time,  the  field  domain 
moves  adiabatically  across  the  whole  sample,  accompanied 
by  the  charge  density  that  fluctuates  in  QWs.  In  this  case,  the 
propagation  of  the  charge-density  fluctuation  pj(t)  —  enj(t) 
can  be  described  by  the  equation2 


(i-pcy-‘-I4D[^(O.^+1(f).4(O]0[4(O]+i>c(i-i>c)/-14?C^o.^i(O.4.(O]«[4(O] 


N 

-pck2l  (i-jPc)*“7'“14d[^(0,^-i(0,4-1(0]^[-4-i(0]-Pc(1-4)^4?[^o.MO,^0] 
xe[-^(0]-/f[My(0,^+1(0,4(0]«[4(0]+^D[40.^-i(0»4-i(0]^[-4->(0].  (3) 


where  the  small  diffusion  current2,8  is  neglected,  0(x)  equals 
one  for  x>0  and  zero  for  j=  1,  2 and  0 ^Pc 
^1  is  the  capture  probability  of  electrons  into  the  QW.  Pc 
depends  on  the  MQW  structure  parameters  and  the  electric 
field  £b  across  the  MQW,  and  can  be  calculated  through  Pc 
=  [lH-Bcoexp(— £cv/£b)]~\  where  and  ^  are  the  high- 
field  capture  coefficient  and  effective  well  capturing  field 
determined  by  the  MQW  structure.1.  In  this  article,  we  will 
simply  take  Pc= 0.5  as  an  example  since  it  does  not  change 
the  qualitative  features  predicted  by  our  model.  The  first  two 
terms  in  Eq.  (3)  represent  the  forward  contributions  from 
capture  current  into  the  jth  QW,  while  the  third  and  fourth 
terms  represent  the  backward  contributions  from  capture  cur¬ 


rent  into  that  QW.  The  last  two  terms  correspond  to  the  for¬ 
ward  and  backward  tunneling  currents  flowing  out  of  the  jth 
QW. 

If  electrons  in  QWs  can  be  approximately  viewed  as  a 
distribution  of  sheet  charges  (zero-thickness),  we  get  the  fol¬ 
lowing  boundary  conditions2  from  the  Maxwell  equations  for 
two  local  fields  on  both  sides  of  the  /th  QW: 

£j(t)-£j-l(t)=  (4) 

e0er 

where  quantum  properties  of  an  electron  gas  in  a  QW  have 
been  ignored,  n2 d  is  the  electron  number  density  in  equilib¬ 
rium,  7=1,  2,..., N,  and  er  is  the  relative  dielectric  constant 
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of  the  well  material  Finally,  the  sum  of  individual  woltage 
drops  on  each  period  (well  plus  barrier)  is  fixed  by  the  volt¬ 
age  Vb(t).  This  restraint  gives  rise  to 
N 

£*5>(0 +(£*+£*) 2  £k(t)=vb{t),  (5) 

k=  1 

and  Lt=LB+N(LB+Lw)  is  the  total  length  of  the  MQW 
structure  in  Fig.  L 

Combining  Eqs.  (3)  and  (4)  for  forward  contributions  we 
know  that  the  sum  of  the  displacement  and  conduction  cur¬ 
rents  is  a  constant  for  Pc=  1,  i.e., 

dS}(t)  9n 

However,  the  conduction  current  itself  is  not  a  constant, 
which  creates  classical  charge-density  fluctuations  and  field 
domains  in  MQWs.  For  non-steady  state,  the  initial  condition 
for  Eq.  (3)  can  be  set  as  p7(0)  =  en2D  if  the  ac  electric  field 
is  applied  to  the  sample  after  f=0.  Equations  (3) -(5)  to¬ 
gether  (totally  2N+ 1  equations)  allow  us  to  simultaneously 
solve  for  the  charge-density  distributions  pk(t)  [or  chemical- 
potential  distributions  fik(t)]  for  k~  1,  2,.. .,7V,  as  well  as  for 
the  local-field  distributions  8k{t)  for  k~0,  1,  2 beyond 
steady  state.  For  steady  state  with  dpj(t)ldt=  0,  only  for¬ 
ward  contributions  will  stay.  In  this  case,  we  only  need  to 
replace  Eq.  (3)  at  each  moment  by5 


D.  Huang  and  D.  A.  Cardimona 

jf[pLj,pLJ+l,£j]  =iVLjW/M  Al 

with y  =  l,  2 Under  this  condition,  the  measured  tun¬ 
neling  current  density  simply  equals  o  >  A  l  >£<)]•  But 

this  does  not  imply  a  uniform  distribution  of  charge  densities 
and  electric  fields.  . 


B.  Quantum  charge  fluctuations 

It  has  been  known  for  a  long  time  that  a  uniform-field 
distribution  will  underestimate  or  overestimate  the 
sequential-tunneling  current  in  MQWs  at  high  or  low  volt¬ 
ages,  respectively,5  while  the  charge  density  will  remain  in 
its  equilibrium  value  even  when  the  uniform  electric  field  is 
an  ac  field.  In  the  classical  field-domain  model,2,5  quantum 
properties  of  electrons  in  QWs  have  been  ignored.  On  the 
other  hand,  it  has  also  been  shown  that  the  quantum  nona- 
diabatic  effects  give  rise  to  a  residual  current  at  zero  ac  ap¬ 
plied  voltage  due  to  a  space-charge-field  effect  in  the  pres¬ 
ence  of  a  uniform  ac  electric  field.3,4,9,10  Therefore  it  is  very 
important  to  include,  simultaneously,  both  the  classical  field- 
domain  effect  and  the  quantum  nonadiabatic  effect  on  the 
sequential  tunneling  of  electrons  in  MQWs  when  an  ac  elec¬ 
tric  field  is  applied  to  the  sample. 

When  the  quantum  nonadiabatic  effect  is  included,  Eq. 
(3)  for  j=  1,  should  be  modified  to 


^^=pc(i-pcy-l^[^o,Mt)A(tms0(t)]+pc^i 

N 

k=j+\ 

lLa\  d£j\t) 


where  the  last  term  represents  the  nonadiabatic  increase  of 
charge  density  in  the  yth  QW,  8^(t)  is  the  nonadiabatic 
space-charge  field  in  the  j th  QW,  S  is  the  cross-sectional 
area  of  the  sample,  and  CgW=(m*e2<S/7r^2) /0[E0 
-/ie(w2D}  T)~\  is  the  quantum  capacitance  with  chemical 
potential  /ie(n2D>  T)  f°r  an  equilibrium  two-dimensional 
electron  gas  in  QWs.  The  nonadiabatic  space-charge  fields 
8^\t)  for  /=1,  2,...,W  in  Eq.  (6)  are  determined  by  the 
following  differential  equations:4,9,10 


d8f{t)  ^  d8j(t) 
Cqw  dt  -CQW  dt 


which  contains  a  “quantum  displacement”  current  due  to 
CQW  as  a  source  term.  The  nonadiabatic  change  of  current 
density  in  Eq.  (7)  is 

A  Jf{t)  ; 

0[  -  £j- j(0]  (8) 
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with  £™+i(t)  =  0.  The  charge-density  fluctuations  in  Eq.  (6) 
now  contain  both  adiabatic  and  nonadiabatic  contributions 
given  by  the  last  term  in  Eq.  (6).  The  total  field  £k(t)  and  the 
average  field  Ej(i )  in  Eqs.  (6)  and  (7)  are  defined  by 


£’k{t) 

£k(t)+£T(t) 

.m. 

.[£*«+ 4-i «]/: 2 

Simultaneously,  Eq.  (4)  should  also  be  modified  to 

(9) 

where  C0—  e0€rS/LB  is  the  classical  geometric  capacitance. 
For  non-steady  state,  the  total  measured  current  density  for 
is  given  by 

^(+)(?)=(i-Pc)^[^o^i(O;£oW]0[foW] 

+  2  {l-Pcf-* 

k=  1 


JLbC0\  d£N(t) 

+[~Y~j  ~dT>  (10) 

which  includes  both  conduction  and  dielectric  displacement 
currents,  where  the  last  term  represents  the  displacement  cur¬ 
rent  from  the  geometric  capacitance.  In  this  case,  the  in¬ 
flowing  injection  current  is  not  equal  to  the  out-flowing  con¬ 
duction  current  from  the  sample  due  to  dpj(t)!dt  #  0.  On  the 
other  hand,  for  Vb(t)<  0  we  have 

+2  (i-^c)*"1 

k=  1 


XJlD[fit(t)-eLB£T(t),fik- 

xei-£Lm+J™[vo,vi(t),£0(tW[£0(t)] 


EftC o\  dS0(i) 
S  j  dt 


(11) 


In  steady  state,  however,  there  is  no  nonadiabatic  charge- 
density  fluctuation,  and  the  total  measured  current  density  is 
mainly  determined  by  the  injection  current,  which  is  given 
for  Vb(t)& 0  by 


which  is  time-dependent  and  different  from  both  C0  and 

^QW- 

From  Eq.  (9)  we  know  that  both  the  field-domain  and 
nonadiabatic  effects  will  cause  charge  density  fluctuations  in 
QWs  in  the  presence  of  an  ac  electric  field.  The  quantum 
capacitance  CQW  only  enters  into  Eq.  (6)  for  the  charge- 
density  fluctuations  but  not  into  Eq.  (10)  for  the  total  mea¬ 
sured  current  density  J*s(±)(f).  Instead,  the  conduction  cur¬ 
rent  is  modified  by  the  nonadiabatic  space-charge  field  £m(t) 
which  is  induced  by  the  “quantum  displacement”  current  as 
shown  in  Eq.  (7).  On  the  other  hand,  the  geometric  capaci¬ 
tance  C0  directly  modifies  the  total  measured  current  in  Eq. 
(10)  as  a  contribution  from  the  dielectric  displacement  cur¬ 
rent  but  does  not  enter  into  the  charge  fluctuations  in  Eq.  (6). 
When  only  the  forward  contributions  are  included,  the  sum 
of  dielectric  displacement  and  conduction  currents  flowing 
into  and  out  of  a  QW  is  a  constant  for  Pc=  1.  Because  the 
conduction  current  flowing  into  the  first  QW  is  simply  the 
injection  current  from  the  contact  layer,  the  change  in  the 
conduction  currents  flowing  into  different  QWs  is  deter¬ 
mined  by  the  variation  of  the  dielectric  displacement  currents 
due  to  the  nonuniform  electric-field  distribution  inside  the 
whole  system.  From  Eq.  (7)  we  further  find  that  even  under 
a  uniform  ac  electric  field,  the  nonadiabatic  conduction  cur¬ 
rent  density  flowing  through  each  QW  is  not  equal  to  the 
sum  of  adiabatic  sequential-tunneling  current  flowing  out  of 
the  QW  and  the  “quantum  displacement”  current  density 
LBCqwd£j(t)/dt  because  d£^{t)f  dt-£  0.  If  the  nonadiabatic 
effect  [or  £”a(f)]  is  small  and  the  backward  contributions  are 
neglected,  Eq.  (6)  can  be  expanded  to  leading  order  as 

dpj(t) 

dt 

j- 1 

~PcS  (i-pcy"*_l4D[p*(o,p*+i(o,4(f)] 


+ pc(  1  -pcy- 1  J^[p0  ,P!  (0  ,£o(0] 


I  s  J  CQW 


dSj(t) 


dt 

?na 

'k  1 


7-1 


+  4  Pc2 


*=i 


£f{t) 


where  the  differential  resistance  in  Eq.  (14)  is  defined  by 


(14) 


J^)=J^o,/*i(0,4(0]  + 


d£dt) 

dt  ’ 


_ 1 


which  can  be  modified  by  the  local  field  at  the  emitter  bar¬ 
rier.  Moreover,  due  to  dpj(t)ldti=  0  in  non-steady  state  we 
can  define  a  total  effective  differential  capacitance  for  the 
MQW  structure 


N 

Oi(/)=2  ci(t)= 

j= i 


dVb(t) 

dt 


-i  at 


<s2 

/=! 


dpjit) 

dt  ’ 


(13) 


The  first  three  terms  in  Eq.  (14)  represent  the  contributions 
from  field  domains.  The  last  two  terms  come  from  the  main 
quantum  nonadiabatic  contributions.  The  other  small  nona¬ 
diabatic  terms  associated  with  the  relative  change  in  charge 
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FIG.  2.  Comparison  of  results  for  total  tunneling  current  /tot(f)  in  a  logarithm  scale  as  a  function  of  the  applied  sinusoidal  bias  from  theory  [in  (a)]  and 
experiment  [in  (b)]  using  ^tot(0 ^ «S[*/2D( 0  4-  AJ^f)].  Here,  barrier  material  is  Al^GaojAs  and  well  material  is  GaAs,  Lw—  50  A,  N=  50,  and  T=40  K.  For 
numerical  calculation,  we  set  nc=6X1017  cm-3,  n2D=4.1X10u  cm""2,  ^=400  s,  <5=7.1  X10-4  cm“2,  and  £0=18  kV/cm.  The  parameters  used  in  the 
experiment  are  nc=lX1018  cm-3,  n2D=2.5X10H  cm“2,  ^=40  s,  S=  1 X  10~4  cm-2,  F0=£0Z,=  6  V,  and  V0  —  £QLt. 


densities  between  neighboring  QWs  are  neglected  in  Eq. 
(14).  Similarly,  Eq.  (7)  can  also  be  expanded  to  leading  order 
as 

dSf(t)  dZj(t)  (  Sf(t)  \ 

dt  dt 

(15) 

which  reduces  to  the  previous  result4,9,10  for  the  case  of  a 
uniform  ac  electric  field. 

111.  NUMERICAL  RESULTS  AND  DISCUSSIONS 

In  this  section,  we  present  numerical  results  for  distribu¬ 
tions  of  both  local  fields  and  charge  densities  in  a  MQW 
sample.  We  first  concentrate  on  effects  of  classical  charge 
fluctuations  described  by  Eqs.  (3)-(5)  for  non-steady  state. 
Alter  this,  effects  of  quantum  charge  fluctuations  will  be  ad¬ 
dressed  based  on  Eqs.  (5)— (9). 

The  sample  we  consider  in  this  article  is  an  AlGaAs/ 
GaAs  MQW  structure.  The  total  number  of  QWs  is  N~  10, 
with  11  barriers.  The  parameters  for  this  sample  are:  well 
thickness  Lw—15  A,  barrier  thickness  LB~ 339  A,  barrier 
height  VQ= 224.5  MeV,  electron  effective  mass  m* 
=  0.065  m0  with  free  electron  mass  ?n0,  electron  areal  den¬ 
sity  n2j)—5X  1011  cm"2,  contact-layer  electron  concentra¬ 
tion  nc=6X  1017  cm-3,  cross-sectional  area  S=  10~4  cm2, 
capture  probability  Pc=0.5,  saturation  velocity  vs=2 
Xi06  cm/s,  saturation  field  £s= 2  kV/cm,  and  relative  di¬ 
electric  constant  €r=  12.  The  ground-state  subband  edge  is 


calculated  to  be  £0  =  44.1  MeV.  The  voltage  is  defined  to  be 
F*(0  =  £c(0£<  with  SM^£0sm(2'TTtftp)  for  t^0,  where 
the  field  amplitude  £0  =  5  kV/cm  and  the  time  period  tp 
=  0.1  s.  The  changes  in  the  sample  parameters  for  the  nu¬ 
merical  calculations  will  be  indicated  in  the  figure  captions. 
[For  typical  MQW  photodetectors,  the  ground  and  excited 
states  are  designed  to  be  bound  and  quasi-bound,  respec¬ 
tively.  However,  for  V0— 224.5  meV  the  choice  of  Lw 
—  15  A  and  LB— 339  A  leads  to  two  lowest  bound  states  in 
QW.  Since  we  only  study  in  this  article  the  effects  of  classi¬ 
cal  and  quantum  charge  fluctuations  on  the  dark  current,  the 
selection  of  excited  state  does  not  change  the  qualitative  fea¬ 
tures  predicted  by  our  model.] 

We  first  show  in  Fig.  2(a)  a  comparison  between  the 
experimental  data  and  the  theoretical  result  from  our  gener¬ 
alized  model  in  Eqs.  (5)-(9)  at  7T=40  K.  This  comparison 
highlights  the  importance  and  uniqueness  of  the  nonadiabatic 
tunneling  effect  included  in  our  generalized  model.  From 
Fig.  2(b)  we  clearly  observe  a  nonzero  residual  current  from 
the  experiment  when  the  ac  voltage  is  swept  through  zero. 
This  phenomenon  cannot  be  understood  at  all  by  Levine’s 
adiabatic  tunneling  model1  which  has  been  widely  applied  to 
study  dark  current  in  quantum- well  photodetectors.  How¬ 
ever,  the  nonzero  residual  current  observed  experimentally 
has  been  -successfully  reproduced  and  explained  by  £™(t) 
#0  due  to  the  nonadiabatic  tunneling  effect  which  is  the 
spirit  of  our  generalized  model  presented  in  Eqs.  (5) -(9)  of 
this  article. 


oJ- 


Barrier  Index  j 


Well  Index  j 


FIG,  3.  Calculated  local  fields  £j(t) 
-£ac(f)  in  (a)  in  different  barrier  lay¬ 
ers  and  density  fluctuations 
—  n2v  in  (b)  inside  different  QWs  at 
times  //^=0.05,  0.11,  0.20,  .and 
0.25  s. 
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.Barrier  Index  j  Barrier  |ndex  j 


FIG.  4.  Calculated  local  fields  Sj{f) 
-£ac(t)  in  different  barrier  layers  at 
the  time  7/^=  0.25  s.  Here,  the  com¬ 
parison  of  £j(t)-EJt)  with  different 
values  of  £0 ,  tp ,  and  T  are  presented 
in  (a)  and  the  comparison  of  those 
with  different  values  of  N,  n2D ,  and 
nc  are  shown  in  (b).  The  changed  pa¬ 
rameters  are  indicated.  The  other  pa¬ 
rameters  used  in  calculations  are  the 
same  as  those  in  Fig.  3. 


A.  Effects  of  classical  charge  fluctuations 

We  show  in  Fig.  3  the  local  fields  £y(f)~£ac(0  in  dif¬ 
ferent  layers  [in  (a)]  and  the  density  fluctuation  rij(t)  -  n1V)  in 
different  QWs  [in  (b)]  at  several  times  t(tp  for  T—ll  K. 
From  (a)  we  find  that  the  field-domain  effect  is  negligible  at 
titp= 0.05  (very  small  applied  field)  because  both  the  injec¬ 
tion  current  and  the  sequential-tunneling  current  jf*  are 
both  extremely  small  in  this  case.  With  the  increase  of  SJj), 
i.e.,  t!tp  increases  from  0.05  to  0.25,  fields  close  to  the  emit¬ 
ter  barrier  are  enhanced  dramatically  relative  to  the  uniform 
field  £ac(0-  This  is  a  result  of  the  huge  current  imbalance 
under  the  uniform  field  £*(0*  as  explained  in  Fig. 
1(b).  At  the  same  time,  fields  close  to  the  receiver  barrier  are 
suppressed  almost  to  zero.  From  Fig.  3(b)  we  find  that  when 
£ac(f)  is  large,  densities  close  to  the  emitter  barrier  are 
greatly  reduced  with  respect  to  the  equilibrium  value  w2D. 
This  is  accompanied  by  a  great  enhancement  of  the  local 
field  at  the  emitter  barrier,  as  shown  in  Fig. 

3(a).  However,  densities  close  to  the  receiver  barrier  remain 
near  to  n2D  due  to  the  suppressed  local  field  £#(0^0-  It  is 
obvious  from  Fig.  3(b)  that  some  electrons  have  been  re¬ 
moved  from  the  sample  since  ^j[nj(t)  —  n2D'i<0,  which  is 
true  even  for  steady  state.  The  calculation  done  here  corre¬ 
sponds  to  a  non-steady  state.  Therefore  the  net  number  of 
electrons  removed  from  the  sample  changes  with  time  under 
an  ac  voltage,  leading  to  a  differential  capacitance  C^(t) 
[see  Eq.  (13)]. 

Figure  4  compares  local  fields  £/(0~£ac(0  at  t!tp 
=0.25  as  a  function  of  barrier  index  j  for  different  values  of 
£o,  i  Tin  (a)  and  different  values  of N,  «2D,  nc  in  (b). 
From  Fig.  4(a)  we  find  that  the  field-domain  effect  is  negli¬ 


gible  at  T'S:  65  K  due  to  very  small  injection  and  sequential- 
tunneling  currents  at  these  temperatures.  The  bigger  the  field 
amplitude  £0  is,  the  larger  £0(t)  will  be.  A  smaller  tp  leads  to 
a  negative  £N(t)  on  the  receiver  barrier  due  to  the  strong 
nonsteady  effect.  This  is  completely  different  from  the 
steady-state  results  in  Refs.  1  and  5  in  which  S^t)  will 
always  be  positive.  Furthermore,  we  find  from  Fig.  4(b)  that 
the  smaller  the  number  of  QWs  N  is,  the  lower  S0(t)  is.  The 
increase  of  n2D  causes  a  larger  sequential-tunneling  current, 
leading  to  a  larger  value  of  £0 (t)  due  to  an  enhanced  current 
imbalance  between  Jfj?  and  .  Conversely,  the  increase  of 
nc  introduces  a  bigger,  injection  current,  leading  to  a  smaller 
value  of  £0(t)  due  to  a  suppressed  current  imbalance  be¬ 
tween  J3^  and  J[D . 

Figure  5  presents  density  fluctuations  ny(r)-n2 d  as  a 
function  of  well  index  j  at  the  moment  of  maximum  £ac(f) 
for  different  values  of  £0 ,  tp ,  T  in  (a)  and  different  values  of 
N,  n2D,  nc  in  (b).  Corresponding  to  Fig.  4(a),  Fig.  5(a) 
shows  a  stronger  charge  depletion  in  the  first  QW  when  tem¬ 
perature  T  is  higher  or  field  amplitude  £q  is  larger.  Accom¬ 
panying  Fig.  4(b),  Fig.  5(b)  also  exhibits  a  larger  deviation  of 
«i(0.  away  from  n2D  for  a  smaller  value  of  nc  or  bigger 
values  of  N  and  n2D . 

As  we  know  from  Eq.  (5),  the  field-domain  distribution 
depends  on  Vb(t)=£ac(t)L,  at  each  moment  t.  Different 
times  imply  different  values  of  local  field  £j(t)  within  the  y'th 
barrier  layer  or  different  values  of  density  rij(t)  inside  the 7th 
QW.  However,  dn dt--  0  for  each  moment  in  steady 
state.  In  non-steady  state,  we  can  define  a  differential  capaci¬ 
tance  as  in  Eq.  (13),  which  is  proportional  to  dn}{t)ldt+  0. 
The  calculated  results  for  Cdi(f)  are  shown  in  Fig.  6(b)  for 


Well  lndex  j  Wall  Index  j 


FIG.  5.  Calculated  density  fluctuations 
nj(t)  ~  n2n  inside  different  QWs  at  the 
time  tltp=  0.25  s.  Here  the  compari¬ 
son  of  rij(t)  —  n1T)  with  different  val¬ 
ues  of  £0 ,  tp,  and  T  are  presented  in 
(a)  and  the  comparison  of  those  with 
different  values  of  N,  n2D,  and  nc  are 
shown  in  (b).  The  changed  parameters 
are  indicated  in  figures.  The  other  pa¬ 
rameters  used  in  calculations  are  the 
same  as  those  in  Fig.  3. 
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FIG.  6.  Effective  total  differential  ca¬ 
pacitance  Cdi(f)  in  (a)  and  compari¬ 
sons  of  conduction  currents  in  (b)  for 
non-steady  state  with  nonuniform  and 
uniform  electric,  fields  and  different 
values  of  £0,  tp,  and  T.  The  conduc¬ 
tion  current  is  calculated  by  setting 
£ja(t)  =  0,  and  the  adiabatic  C^r)  is 
also  calculated  with  £“(?)  =  0.  The 
changed  parameters  are  indicated  in 
figures.  The  other  parameters  used  in 
calculations  are  the  same  as  those  in 
Fig.  3. 


different  values  of  £0,  tp,  and  T.  It  is  clear  from  (b)  that 
Cdi(0  is  always  negative,  increases  greatly  with  £ac(0>  and 
becomes  nearly  zero  at  T^65  K.  The  value  of  Cdi(f)  is 
dominated  by  Cdl(t)  in  the  first  QW.  Whenever  \nx(i) 

—  «2dI  is  larger  in  Fig.  5(a),  —  Cdi(/)  becomes  bigger  in  Fig. 
6(b).  From  an  experimental  point  of  view,  the  conduction 
current  detected  by  the  receiver,  as  described  by  Eq.  (14) 
without  the  last  term,  is  of  great  interest.  We  show  in  Fig. 
6(a)  the  conduction  current  as  a  function  of  time  tltp .  In 
order  to  see  the  effects  of  classical  charge  fluctuations,  we 
compare  the  conduction  current  calculated  from  Eq.  (14)  by 
excluding  the  last  term  and  setting  £™(t)  =  0  (solid  curve) 
with  the  current  calculated  from  Eq.  (2)  by  setting  £k(t) 

—  £ac(0  for  a  uniform  electric  field  and  nk(t)  —  n2x>  (dashed 
curve).  When  £ac(f)  is  small  with  t/tp^0.05y  solid  and 
dashed  curves  agree  very  well  due  to  £jv(f  )**£«?( 0  as  seen 
from  Fig.  3(a).  However,  the  solid  curve  falls  below  the 
dashed  curve  after  tltp  — 0.05  and  then  changes  sign  when 
t/tp> 0.1.  The  solid  curve  in  (a)  can  qualitatively  be  ex¬ 
plained  by  Fig.  3(a),  where  with  the  increase  of  tltp ,  £^{t) 
starts  at  £ac(f),  and  then  becomes  negative,  and  finally  ap¬ 
proaches  zero.  The  negative  conduction  current  observed  in 
(a)  is  a  direct  result  of  our  non-steady-state  conditions,  which 
is  completely  different  from  those  in  Refs.  1  and  5  for  steady 
state. 

To  summarize  the  above  observations,  effects  of  classi¬ 
cal  charge  fluctuations  only  play  a  role  under  the  condition  of 
a  large  imbalance  between  the  injection  and  sequential- 
tunneling  currents  near  the  emitter  barrier.  This  can  be 
achieved  by  increasing  either  the  temperature  T  or  the  field 
amplitude  £0  * 

B.  Effects  of  quantum  charge  fluctuations 

Effects  of  classical  charge  fluctuations  depend  on  the 
geometric  capacitance  C0,  as  shown  by  Eq.  (4).  However, 
the  quantum  capacitance  CQW  starts  to  play  a  role  when  the 
nonadiabatic  effect  is  included,  as  seen  from  Eq.  (9).  The 
striking  thing  is  that  the  effects  of  quantum  charge  fluctua¬ 
tions  become  more  and  more  important  as  T <  65  K,  while 
the  effects  of  classical  charge  fluctuations  are  negligible  at 
these  temperatures.  We  compare  in  Fig.  7  two  calculated 
nonadiabatic  fields  £™(t)  at  tltp=  0.25  as  a  function  of  well 
index  j  for  T=40  and  77  K.  From  the  figure  we  know  that 
£™(t)  decreases  with  j  due  to  field-domain  effects  at  T 


=77  K.  However,  £j\t)  becomes  independent  of  j  at  T 
=  40  K.  More  importantly,  £]\t)  increases  with  reducing  T. 

The  effect  of  nonadiabatic  field  £Ja(0  can  be  seen  more 
clearly  from  the  calculated  local  fields  £j(t)~ £j^t)  in  dif¬ 
ferent  barrier  layers  in  Fig,  8(a)  and  density  fluctuations 
nj(t)  —  n2 d  in  different  QWs  in  Fig.  8(b)  at  */rp= 0.25  with 
various  values  of  £0,  tp ,  and  T .  By  comparing  solid  and 
dashed  curves  in  (a)  we  find  that  £0(t)  is  reduced  by  a  factor 
of  2  as  £0=5kV/cm,  ^=0.1  s,  and  T=77K  (squares). 
Moreover,  £0(t)  decreases  even  more  due  to  £]\t)  when  tp 
is  reduced  to  0.05  s  (stars)  due  to  stronger  nonsteady  effects, 
but  it  decreases  much  less  when  T  is  reduced  to  40  K  (tri¬ 
angles)  due  to  smaller  injection  and  sequential-tunneling  cur¬ 
rents.  However,  £0(t)  is  enhanced  by  nonadiabatic  effects 
when  £0  is  reduced  to  1  kV/cm  (circles).  All  of  these  features 
have  a  strong  influence  on  C^(f),  as  discussed  in  Fig.  6(a). 
On  the  other  hand,  from  (a)  we  also  find  that  £N{t)  becomes 
a  much  more  negative  value  for  £0  =  5  kV/cm,  tp— 0.1  s,  and 
T—ll  K  when  £™{i)  is  included  in  our  calculations  (solid 
curve  with  squares)  than  when  it  is  excluded  (dashed  curve 
with  squares).  The  same  situation  occurs  when  either  tp  is 
reduced  to  0.05  s  (stars)  or  £0  is  reduced  to  1  kV/cm  (circles) 
but  with  a  smaller  overall  magnitude  compared  to  the  curves 
with  the  squares.  As  explained  in  Fig.  6(b),  the  features  ob¬ 
served  for  £x(t)  in  Fig,  8(a)  will  affect  the  conduction  cur- 


FIG.  7.  Calculated  nonadiabatic  space-charge  fields  £]\t)  at  tftp-  0.25  as  a 
function  of  well  index  j  for  T~  40  and  77  K  (left  scaled),  respectively.  The 
other  parameters  used  in  calculations  are  the  same  as  those  in  Fig.  3. 
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FIG.  8.  Calculated  local  fields  Sj(t) 
-£ac(f)  in  different  barrier  layers  and 
density  fluctuations  7ij(t)-n2 D  inside 
different  QWs  at  tltp= 0.25  with  dif¬ 
ferent  values  of  £0  (black/red  curves), 
tp  (black/blue  curves),  and  T  (black/ 
green  curves).  are  com¬ 

pared  in  (a)  'with  (solid  curves)  and 
without  (dashed  curves)  £“(/),  and 
the  similar  comparisons  of  ny-(f) 
~~n2D  are  presented  in  (b).  The 
changed  parameters  are  indicated  in 
figures.  The  other  parameters  used  in 
calculations  are  the  same  as  those  in 
Fig.  3. 


rent  detected  at  the  receiver  layer.  Effects  of  quantum  charge 
fluctuations  are  reflected  in  the  calculated  n;-(f)  -  n2v  in  (b), 
where  the  reduction  of  charge  in  the  QWs  is  greatly  in¬ 
creased  except  for  the  case  with  a  small  value  of  £0  (circles). 
Although  there  is  a  strong  dependence  of  £j(t)  on  index  j 
near  the  emitter  barrier,  there  is  very  little  dependence  of 
rijit)  on  j  there  after  £]\t)  is  included  in  the  calculations. 

IV.  CONCLUSIONS 

In  conclusion,  by  including  nonadiabatic  space-charge- 
field  effects  we  have  generalized  the  previous  theories  for 
studying  field-domain  effects  in  MQW  photodetectors  in  the 
presence  of  an  ac  voltage.  We  have  found  from  our  numeri¬ 
cal  calculations  that  field-domain  effects  are  only  important 
at  high  temperatures  or  high  voltages,  which  implies  the  ex¬ 
istence  of  significant  injection  and  sequential-tunneling  cur¬ 
rents  in  the  system.  We  have  further  found  that  nonadiabatic 
effects  become  much  more  visible  at  low  temperatures  and 
low  voltages  when  the  field-domain  effects  are  negligible. 
The  time  duration  for  nonadiabatic  effects  has  been  found  to 
relate  to  the  quantum  capacitance,  while  the  dielectric  dis¬ 
placement  current  has  been  shown  to  be  controlled  by  the 
geometrical  capacitance.  For  non-steady  state,  we  have 
found  a  negative  differential  capacitance  in  the  system,  and 
we  have  predicted  a  negative  conduction  current  under  a 
positive  voltage.  f 

In  this  article,  we  have  assumed  that  the  capture  prob¬ 
ability  is  independent  of  electric  field.  This  can  be  justified 
by  the  fact  that  the  capture  probability  is  near  a  constant  at 
low  electric  fields.  In  the  presence  of  incident  photons,  the 
conduction  current  flowing  through  the  MQW  sample  will  be 
the  sum  of  sequential-tunneling  and  photoexcited  currents. 
From  our  studies  in  this  article,  we  predict  that  the  field 
domain  effects  which  are  significant  at  high  temperatures  or 
high  photon  fluxes  will  strongly  affect  both  tunneling-  and 
photo-currents.  On  the  other  hand,  nonadiabatic  effects 
modify  these  currents  at  low  temperatures  and  low  photon 


fluxes  when  the  current  flowing  through  MQWs  is  large. 
These  latter  conditions  are  of  utmost  importance  in  any 
space-based  detector  applications. 

One  of  the  main  results  of  this  work  is  the  prediction  of 
the  dominance  of  field-domain  effects  at  77  K  and  above  and 
the  dominance  of  nonadiabatic  effects  at  40  K  and  below. 
The  experimental  observation  of  anomalous  nonadiabatic  be¬ 
havior  in  sequential-tunneling  current  at  low  temperatures 
and  the  suppression  of  it  at  high  temperatures  has  already 
been  demonstrated.2’3  For  a  bound-to-bound  multiple- 
quantum-well  photodetector,  if  there  exists  only  one 
negative-differential-conductance  peak  at  a  voltage  predicted 
by  a  uniform-field  model  at  low  temperatures,  it  is  an  indi¬ 
cation  of  the  suppression  of  the  field-domain  effect.11 
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